Bilayer manganites: polarons in the midst of a metallic breakdown 
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The exact nature of the low temperature elec- 
tronic phase of the manganite materials family, 
and hence the origin of their colossal magnetore- 
sistant (CMR) effect, is still under heavy debate. 
By combining new photoemission and tunneling 
data, we show that in La2-2 X Sri_|_2 X Mn2 07 the po- 
laronic degrees of freedom win out across the 
CMR region of the phase diagram. This means 
that the generic ground state is that of a sys- 
tem in which strong electron-lattice interactions 
result in vanishing coherent quasi particle spec- 
tral weight at the Fermi level for all locations in 
k space. The incoherence of the charge carriers 
offers a unifying explanation for the anomalous 
charge-carrier dynamics seen in transport, op- 
tics and electron spectroscopic data. The stack- 
ing number N is the key factor for true metal- 
lic behavior, as an intergrowth-driven breakdown 
of the polaronic domination to give a metal pos- 
sessing a traditional Fermi surface is seen in the 
bilayer system. 

Competition between local lattice distortions leading 
to anti-ferromagnetic, charge and orbital (CO) ordering 
on the one hand, and mixed valence character promot- 
ing metallic ferromagnetic double exchange on the other, 
determines the transport transport properties [1] of the 
manganite materials family and is proposed to lie at 
the root of their CMR effect [2]. (La,Sr) N+ iMn N 3N +i, 
abbreviated LSMO (where N is the number of stacked 
MnC>2 planes between [La,Sr]0 block layers [3]) displays 
a remarkable decrease in metallic character with decreas- 
ing N [4], see Fig. SI. Metallic behavior thrives in the 
phase diagram of cubic LSMO [5], whereas the bilayer 
analogue is metallic only in a narrow Sr-doping and tem- 
perature regime [6] , giving rise to the largest CMR effect 
[7]. The more strongly 2D, single layer compound shows 
neither metallic nor CMR behavior [8] . 

Focusing on bilayer LSMO within the CMR-region of 
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the phase diagram, the prevailing picture from structural 
studies is one of polarons existing above To On cool- 
ing towards To these short range versions of the CO 
order typical of the insulating compositions, become in- 
creasingly correlated [9, 10]. Eventually double exchange, 
leading to an itinerant, metallic state, takes over. 

This metallic state for x — 0.4 has been shown to 
support small quasi-particles (QPs) in the spectral func- 
tion measured by angle-resolved photoemission (ARPES) 
[11]. These signal coherent electronic excitations, albeit 
strongly dressed with lattice distortions, and are seen as 
evidence for a novel and elusive state of matter known as 
a polaronic metal [11, 12]. 

Other ARPES studies paint a different picture, with 
stronger QP features observed at low T that persist up to 
temperatures of order 1.5Tc [13-15], despite the system 
being nominally insulating. In contrast, scanning tunnel- 
ing microscopy / spectroscopy (STM/S) studies reported 
gaps in the local density of states near-E F for x = 0.30 
[16] and 0.325 [17], both in the metallic and insulating 
temperature regimes. Finally, new neutron diffraction 
data for x = 0.4 has shown that even far below Tc - at 
10 K in the metallic state - polarons remain as fluctua- 
tions that strongly broaden and soften phonons near the 
wave vectors where the charge order peaks would appear 
in the insulating phase [18]. 

Here, a combination of ARPES and STM/S reveals 
that the bilayer manganites still have a number of tricks 
up their sleeves. Firstly, we show that the intrinsic spec- 
tral response of these systems is pseudogapped, with neg- 
ligible coherent spectral weight at Ep anywhere in k- 
space at any temperature, across the CMR-region of the 
phase diagram. Secondly, we show that the strong QP 
features seen in ARPES studies - also above Tc - are 
due to the unavoidable presence of TV ^ 2 stacking-fault 
intergrowths. 

These new insights clear the way for a unified interpre- 
tation of the physical properties of CMR bilayer LSMO in 
terms of strongly incoherent charge carriers, and suggest 
that local control of the dimensionality of such mangan- 
ites - via tuning the stacking number, N - may offer a 
novel route to new functional nanostructures. 



2 




Figure 1: STM/S on bilayer manganites. (a) Large field of view of a region with (an above average number of) half unit 
cell steps at the surface and (b) corresponding line scan (x = 0.3, V = —500 meV, I — 20 pA). (c) Semi-ordered square-like 
texture for x=0.36, (V = —450 meV, I = 150 pA), which is more clearly seen in the line trace (d) through the center of 
the autocorrelation (inset of d, 220x220A 2 ). A correlation length of ~2.3 nm is determined from fitting Lorentzians to the 
autocorrelation features, (e) STS spectra - I(V) curves - showing a gap of more than 100 meV for all doping concentrations 
studied (V = 500 meV, I ~ 50 pA). The spectra offset with respect to each other for clarity, indicated by the dashed lines, 
(f) Comparison between symmetrised, k-integrated ARPES spectra and STM dl/dV traces, (g) Topograph displaying atomic 
resolution (V = 500 meV, I = 140 pA) with (lxl) and v2 x y2 spots in the Fourier transform (see inset, black and white 
arrows respectively) which are not found in the corresponding low energy electron diffraction images, (h) Topography of stepped 
surface taken on the same region as (g), where the step height corresponds to that expected for single layer LSMO. 



Imaging the cleavage surface of a typical bilayer LSMO 
single crystal in the CMR doping region using low tem- 
perature STM/S yields large, flat terraces, see Figs, la 
and b. One surprising characteristic of these STM data is 
that the surface atoms have proven extremely difficult to 
image. In the literature, the only reported observations of 
atomic resolution in bilayer LSMO concerned small nm- 
sized patches [16, 17]. The terraced, flat and debris-free 
surfaces we image may lack atomic corrugation, however 
they do posses a spatial texture in the tunneling signal. 
These structures are usually disordered, but are - in some 
cases - ordered into a semi-regular, square-like lattice, 
displaying characteristic length scales of order 2-3 nm, as 
can be seen in Fig. lc. The quasi-ordered nature of these 
regions at the surface is evident from the autocorrelation 
traces from the STM topographs, as shown in Fig. Id, 
which display clear structures at distances between 5-15 
units cells. 

Comparison with low energy electron diffraction 
(LEED) data recorded in the STM chamber shows that 
the orientation of the quasi-periodic structures is not 
pinned to the underlying atomic lattice. The tunneling 
spectra of bilayer LSMO from all cleavage surfaces mea- 
sured are gapped (symmetrically) around Ep over an en- 
ergy range of several 100 meV at low temperature for all 



doping concentrations studied, as shown in Fig. le. The 
tunneling spectra themselves show no major variations 
from point to point. This suggests that the observed spa- 
tial textures seen in the STM data are a subtle, 'higher 
order' structural ordering phenomenon, leading to mod- 
est, periodic signatures, somewhat reminiscent of the 2.5 
nm correlated regions seen in the high-T, polaronic liq- 
uid phase described in Ref. [10]. The authors of Ref. [2] 
suggest that aggregates of polarons resembling a charge 
ordered state should be observable in STM experiments 
above Tq . Our new data show that these kind of entities 
could exist at the surface even at low temperatures, pos- 
sibly being a log-jammed form of the fluctuating polarons 
observed in recent neutron experiments [18]. 

During our STM investigation of numerous cleavage 
surfaces of bilayer LSMO, very clear atomic resolution 
was found over a large terrace of one particular cleave, 
(shown in Fig. lg). Analysis of the Fourier transform 
shows reduced in-plane symmetry, with clear y2 x y2 
spots shown in the inset to Fig. lg. This is atypical 
for bilayer LSMO, which has a (1 x 1) tetragonal lattice 
symmetry in the ferromagnetic 'metallic' phase that we 
observe in all our LEED data (see Fig. S2c). An impor- 
tant observation linked to this large scale STM imaging 
with atomic contrast is that the step bordering this ter- 
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race is only one quarter of the N — 2 c-axis unit cell in 
height (5 A, see Fig. lh). Additionally, the STS spectra 
from this region were more strongly gapped compared 
to those from non-atomically resolved regions displaying 
10 A step heights, as shown in Fig. If. Taken together, 
these facts form a compelling argument that in Figs, lg 
and h we are - in fact - imaging a stacking fault in the 
Ruddlesden-Popper manganite in which an extra La202 
block (or blocks) creates a region of the bilayer crystal 
which is effectively single layer (La,Sr) 2 Mn04. The fact 
that atomic contrast is readily observed in STM of N = 1 
LSMO [19] also further strengthens our interpretation. 

The observation with STM/S of such a surface inclu- 
sion with a different stacking number N is rare, but is 
not, in itself, wholly surprising. Numerous /iSR studies 
[20], magnetization measurements [21, 22] and transmis- 
sion electron microscopy studies [23, 24] have established 
that stacking faults occur even in the very best crystals at 
the w 1% level [21, 22]. These intergrowths vary - locally 
- the stacking number N and are the cause of anomalous 
steps in the magnetization above the bulk Tc common 
in bilayer LSMO between 200 and 350 K: small patches 
with N values above two deliver a higher Tc- Due to 
the strong connection between magnetism and metallic 
behavior in the manganites, it is a simple step to reason 
that N > 2 intergrowths would also be metallic above 
and beyond the bilayer Tc- In contrast, an TV = 1 in- 
tergrowth will be more insulating than bilayer LSMO, in 
keeping with the STM/S data shown in Fig. If. 

Armed with a heightened awareness that (unavoidable) 
intergrowths in bilayer LSMO are present in real samples, 
and that signals from these regions can be picked up in 
spectroscopic experiments, we now turn to our ARPES 
investigations. In total, we have measured dozens of bi- 
layer LSMO single crystals from two different sources 
with Sr doping levels spanning the entire CMR region 
of the phase diagram. ARPES spectra were taken at 
low temperature over a square grid across the samples in 
steps of 100 fim, see Fig. 2b. Using such a procedure, two 
types of spectral signature are observed on different loca- 
tions: (i) from only 1 — 5% of the cleavage surface, sharp 
QP-peaks are observed at all /c-locations on the Fermi 
surface, as shown in Fig. 2c and (ii) vanishingly small 
spectral weight is observed near Ep from the majority of 
the sample surface, as can be seen in Fig. 2a for 10 K and 
in the comparison made in Fig. 2d. Fig. 2 also shows the 
ARPES signatures from the non-QP-peaked and QP- 
peaked regions for T>Tq. For the non-QP-peaked re- 
gions, a modest shift in spectral weight to higher binding 
energies can be seen. For the minority, QP-peaked loca- 
tions, the near Ep spectral weight is reduced at 150K (a 
temperature that is almost twice the bilayer Tc for this 
doping), however the /(-E)-traces in Fig. 2d still display 
considerable intensity at Ep , and certainly no signs of a 
pseudogap. This despite the fact that the vast majority 
of the sample is in the paramagnetic, insulating regime. 
In Figs. S3 and S4, we show that for varying doping lev- 
els the Fermi surface (FS) topology (i.e. the number of 



bands), the FS area and the QP-intensity vs. tempera- 
ture for the QP-peaked regions are not in line with the 
bulk properties of bilayer LSMO. Additionally, when ob- 
served for a particular doping, the QP-peaks display a 
cleave-to-cleave variation in the number of FS sheets ob- 
served around the X-point ranging from one to four. 

Thus, combining our atomic scale tunneling data with 
these ARPES data it is clear that the strongly QP- 
peaked spectral signature is not representative of bilayer 
LSMO. Instead, the intrinsic signature of this material is 
that of a pseudo-gapped, very bad metal that supports 
vanishingly small QP spectral weight at Ep at any Al- 
location. This may sound surprising for a metal, but the 
bilayer manganites are far from being normal metals. In 
the 'metallic' phase, the ab-plane resistivity is in excess of 
the Mott maximum of 10~ 3 flcm [26, 27], which is quite 
unlike the situation in their cubic cousins (see Fig. SI). 
In addition, the optical conductivity of bilayer LSMO 
shows an incoherent Drude peak, even down to the lowest 
temperatures [28, 29], again, unlike the situation in the 
N = oo, cubic materials [30]. Therefore, it is clear that 
the strongly QP-peaked regions of the cleaves are simply 
signals from stacking faults [22] located sufficiently close 
to the crystal termination (Fig. 2e and 2f) which conse- 
quently contribute strongly to the near-Ep photoemission 
signal. In this manner, the anomalous temperature de- 
pendence of the QP-peaks (Fig. S4c) also falls naturally 
into place. Fig. S5 provides additional arguments involv- 
ing the orthorhombic crystal symmetry of the stacking 
fault regions observed in ARPES data; here we note that 
the inset to Fig. lg already shows the orthorhombic sur- 
face of an intergrowth (in this case N = 1). 

The new experimental data from real space and k— 
space probes presented above provides a unifying frame- 
work in which to understand all the published ARPES 
data, which is a welcome simplification. However, the bi- 
layer manganites still present a richness - in particular in 
the form of a highly sensitive doping dependence - which 
is a major challenge to our theoretical understanding. Al- 
though in our experiments we have adopted a wide range 
of experimental conditions including geometries and po- 
larizations [31] very similar to those of Rcfs. [11] and [12], 
we have not observed the nodal quasiparticles reported 
for the x = 0.4 composition. This may seem puzzling un- 
til one realizes that the phase diagram of bilayer LSMO 
possesses numerous line phases. Various narrow regions 
in the phase diagram show special types of magnetic 
and/or orbital ordering: x — 0.30 (AFM metal), 0.50 
(CE-type charge and orbital ordering) and 0.60 (AFM 
metal) [32-34]. Remarkably, deviations of only 0.01 in 
doping level (e.g. away from x = 0.60 [34, 35]) have a 
major effect on the electronic behavior. Accordingly, we 
suggest that exactly at or very close to a composition 
of a; = 0.40 there is also such a line phase, whose ex- 
act composition has not been 'hit' in the ARPES data of 
Refs. [13-15, 36, 37]. 

Our data - covering samples right across the CMR re- 
gion of the phase diagram - shows that the generic elec- 
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Figure 2: Spatial ARPES mapping of bilayer LSMO. (a) Typical I(E, fc)-images taken around the (7r,0) point of fc-space for 
temperatures well below and above Tc for an x — 0.425 sample from a non-QP-peaked region of the cleave (hi/ = 56 eV). A 
clear absence of QP features characterizes the data both above and below the bulk Tc- (b) Optical micrograph of the post- 
cleavage surface of a bilayer LSMO crystal (from an x — 0.36 sample). At each red grid point within the blue box an ARPES 
spectrum has been taken. The synchrotron spot size is shown to scale. The regions which do support QP-peaked spectra are 
indicated by the green shaded areas. The proportion of the surface showing QP-peaked spectra is unusually abundant for this 
sample. The inset below panel b shows a schematic of the Brillouin zone, after Ref. [25], with the ARPES cuts shown in a and 
c indicated in blue, (c) Analogous data to those of panel a, now for a QP-peaked region of the same sample, (d) /(S)-curves 
taken along the dashed lines in (a) and (c). For the QP-peaked region, there remains significant spectral weight at E_f for 
150K, despite the fact that T — > 2 Tc- In sharp contrast, the non-QP-peaked I{E) curves display a pseudogap-like lack of 
EF-spectral weight both below and above Tc- Panel (e) shows the explanation for the striking lateral position dependence of 
the electronic structure. Depicted is a schematic side-on view of a cleaved bilayer LSMO crystal with a stacking fault (illustrated 
is an N — 5 intergrowth) just under the crystal surface. Such inclusions support QP-states to temperatures well above the 
bilayer Tc- The intrinsic spectral signature of bilayer LSMO is that seen in panel a and the blue traces in panel d. (f) In 
connection to the schematic in e: a transmission electron microscopy image, adapted from [22], showing an TV = 5 inclusion in 
a bilayer LSMO single crystal. 



tronic signature of bilayer LSMO is that of a pseudo- 
gapped, very bad metal with vanishingly small QP spec- 
tral weight at the Fermi level. Hence, 'nodal metallicity', 
although of great interest in the context of comparisons 
with the high T c superconducting cuprates, is apparently 
not a necessary factor for bringing about CMR behavior, 
and appears to be confined to the doping level x = 0.40. 

Before reaching a final conclusion as regards the 
physics behind a metal which has QP spectral weight be- 
low our detection limit, one further issue needs to be dealt 
with, and that is whether the surface of bilayer LSMO 
is even less like a regular metal than the bulk. X-ray 
resonant diffraction has shown the first pair of Mn02~ 
planes nearest to the surface of bilayer LSMO cleaved in 
air [38, 39] or in vacuum [40] to possess no ferromagnetic 
order and thus the surface could be expected to display 



non-metallic behavior. A recent (hard) x-ray photoemis- 
sion study of the same crystals as studied here found no 
major difference between the surface and bulk in terms of 
charge transfer or composition [41], thus the differences 
between the outermost bilayer and the rest are subtle in 
nature. While recent LEED studies [40] have given no 
evidence for a lowering of the 2D surface symmetry, for 
example via a surface reconstruction, a finding our own 
LEED data support, a contraction of the apical bond 
length for the outermost Mn02-plane was observed [40] , 
which could impact the mobility of the charge carriers at 
the surface. 

Here it is essential to realize that the photocurrent 
measured in ARPES certainly originates from deeper in 
the crystal than the first bilayer only, even at excitation 
energies in the vacuum ultraviolet. Therefore, we can 



5 



conclude that the bad-metallic, fully pseudo-gapped be- 
havior we observe is indeed a genuine characteristic of the 
bulk for the vast majority of the CMR-compositions in 
bilayer LSMO. Although disentangling the spectral sig- 
natures of the outermost and deeper lying bilayers is dif- 
ficult, we do note that the energy range at Ep over which 
the spectral weight is suppressed in ARPES is smaller 
than the negative bias gap seen in STM/S, the latter 
probing strictly the electronic states at the surface. This 
is evident from Fig. If where STS spectra from bilayer 
LSMO and from the N = 1 intergrowth are compared 
with symmetrized, k-integrated photoemission spectra 
from ARPES maps (such as those in Fig. S3) for a QP- 
peaked and a non-QP-peaked region of the surface. 

We now return to the key issue of the physical na- 
ture of the charge carriers responsible for electronic trans- 
port of the bilayer manganites. The most suitable pic- 
ture - considering the results on both the dynamics of 
the charge carriers and the subtle yet clear structural 
anomalies they bring with them - is that of polaronic 
charge carriers above Tq, but also in the ferromagnetic 
(bad) metallic state. The fact that the STM signatures 
of the weakly self-organised polaronic carriers when 'log- 
jammed' at the surface resemble the polaronic correla- 
tions seen above Tc suggest that in the bulk of the bilayer 
LSMO systems - which is more metallic than the surface 
- the system teeters on the edge of a breakdown of the 
fragile polaronic metal state into an insulting, charge and 
orbitally ordered state. This borderline situation, cou- 
pled to the disorder and enhanced fluctuations present 
in these quasi-2D systems, delivers all the ingredients for 
the colossal magnetoresistance transition [2] . 

We conclude this paper by making a connection back to 
the relationship between the materials physics and crys- 
tal chemistry of the manganite family. The physics de- 
scribing the manganites depends crucially on the propen- 
sity of the system to form ordered textures in spin, charge 
and orbital occupation. For the polaronic metal state to 
remain stable, it is vital that the degeneracy in the e g 
orbital manifold is preserved, as this is the ticket to the 
double-exchange energy reduction that encourages hop- 
ping of the carriers. This condition - in turn - stipulates 
equality in the equatorial and axial bond lengths of the 
Mn06 octahedra. The deviating axial bond lengths at 
the surface of bilayered LSMO compared to the bulk, 
as observed in Ref. [40] lift this degeneracy and thus 
push the surface of this material further into the insu- 
lating regime than bulk: the correlated polarons can be- 
come static, as we observe with our STM measurements. 
On the other hand, the N > 2 inclusions we observe, 
structurally bear greater resemblance to the cubic com- 
pound, and thus show a higher propensity towards metal- 
lic behavior, including the existence of coherent spectral 
weight. Paradoxically, the fragility of the polaronic metal 
state in the bilayer systems is also the key to their colossal 
magnetoresistance, as it delivers the precarious balance 
between weakly metallic and insulating behavior required 
for such an enormous sensitivity to the extra impulse pro- 



vided by an external magnetic field. 

The results presented here formulate a clear challenge 
to develop a general theory for the transport in N = 2 
systems involving practically incoherent charge carriers - 
fluctuating polarons - while also capturing the sensitiv- 
ity to the stacking number, N. Aside from the lattice 
polaron generally considered, different types of polaron, 
such as spin and orbital polarons, should be taken into 
consideration [42-44]. An interesting proposal for the 
charge carrier dynamics in the manganites is the exis- 
tence of so called Zener polarons [45, 46]. In such a 
polaron, the charge carrier is not localised on a single 
manganese atom, but on two neighbouring manganese 
atoms which are ferromagnetically coupled by the Zener 
double-exchange mechanism. It has been shown that all 
these different types of polaron lead to a large incoherent 
spectral weight such as is seen in our ARPES investiga- 
tions [42-44, 47]. 

The new insight we have gained points towards the 
great potential of e.g. layer-by-layer thin-film engineer- 
ing to generate tailor-made heterostructures, not only 
to lead to enhanced transition temperatures [48], but 
in combination with modern lithographic and patterning 
methods to tune and improve magnetoresistive properties 
on the sub-micron scale in a new generation of complex 
oxide devices. 



I. METHODS. 

STM. The STM data were recorded using a commer- 
cial variable temperature UHV microscope from Createc 
GmbH. W tips were prepared by electrochemical etching 
followed by in-situ conditioning before each measurement 
using a Au(788) single crystal. In all cases, the spectral 
shapes obtained were independent of the tip to sample 
distance, a sign of good vacuum tunneling conditions. 
All investigated samples have been cleaved in situ prior 
to measuring. The crystallographic orientation was de- 
termined directly after the STM measurements using in 
situ low energy electron diffraction. The setup voltages, 
as well as the tunneling currents are indicated in the fig- 
ures. 

ARPES. The angle-resolved photoemission data 
presented here were obtained using the following beam- 
lines and end-stations: i) the UE112-PGMa beamline at 
the Hclmholtz Zentrum Berlin (BESSY II storage ring), 
Berlin, Germany, coupled to an SES100 analyzer; ii) the 
UE112-PGMb beamline at BESSY coupled to the R4000 
analyzer of the leubed end-station and iii) the SIS beam- 
line at the Swiss Light Source (SLS), Villigcn, Switzer- 
land, equipped with an SES2002 analyzer. The total ex- 
perimental energy broadening at 20 K was set to 30 meV, 
20 meV and 15 meV, for the three end-stations, respec- 
tively. The momentum resolution was 0.027r/a along the 
analyzer slit at the excitation energies used. The pho- 
ton energies used arc indicated in the relevant figures. 
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High quality single crystals of LSMO were grown in both 
Oxford and Amsterdam using optical floating zone tech- 
niques, and were cleaved in situ and characterized us- 
ing LEED before ARPES measurements were conducted. 
For further details concerning the sample quality and 
characterization, see Fig. S2. 
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I. DEPENDENCE OF THE TRANSPORT PROPERTIES ON THE DIMENSIONALITY 

Figures Sla-c show the phase diagrams taken from literature of (La,Sr)jv + iMnjv03Ar + i, for N = 1 [1], 2 [2, 3] and oo 
[4] . Comparison of the diagrams clearly illustrates the strong dependence of the degree of apparent metallic character 
on the number of stacked Mn02 planes, N. Whereas the cubic compound is truly metallic over a wide range of Sr 
doping content, the N = 1 composition is not metallic for any doping concentration. The N = 2 bilayer material is 
somewhere in between these two extremes, displaying metallic transport, but with an absolute value of resistivity in 
excess of the Mott maximum. To stress this latter point, Fig. Sib and Fig. Sic show the temperature dependence of 
resistivity for various doping levels of the cubic compound [5] and the x = 0.4 doping concentration of the bilayered 
LSMO material [6], respectively. As can be seen, the absolute value of resistivity of the cubic composition on the 
border of the insulating phase (x — 0.175) is still considerably less than the 'metallic' x = 0.4 bilayered material. 
This is also apparent from optical conductivity measurements, where Drude weight is absent in the bilayered material 
[7, 8] , while it is clearly present in the cubic compound [9] . 




Doping (x) Doping (x) Doping (x) Temperature (K) 

S 1: Magnetic and transport properties of the LSMO Ruddlesden-Popper series. (a) Magnetic phase diagram of 
(La,Sr)jv+iMnjv03Ar+i as a function of Sr doping for the perovskite analogues N = 1 (left, after Ref. [1]), the bi-layered 
N = 2 (middle, after Ref. [2, 3]) and single layered N = 1 (right, after Ref. [4]). (b) Resistivity as function of temperature 
for various metallic compositions of perovskite LSMO (N = oo, after Ref. [5]). (c) Resistivity as function of temperature for 
bi-layered LSMO with x = 0.4 (from Ref. [6]). 



2 



II. SAMPLE AND CLEAVAGE SURFACE QUALITY OF THE BILAYER LSMO CRYSTALS 

The samples measured throughout this study were grown using the optical floating zone technique in both Oxford 
(x = 0.30, 0.325, 0.35, 0.375, 0.40, 0.425, 0.475) [10, 11] and Amsterdam (x = 0.30, 0.36, 0.40). Electron probe 
micro-analysis measurements on selected LSMO samples revealed that the composition of the single crystals was 
within a few percent of the nominal doping concentration and did not vary substantially across the crystal surface. 
All measured crystals were first characterised by magnetometry to ensure the sample has a sharp magnetic transi- 
tion. Figure S2a shows a typical magnetization versus temperature trace, showing a sharp transition into the low 
temperature ferromagnetic phase. Figure S2b shows the phase diagram constructed from the transition temperatures 
of all doping concentrations used in this study, plotted on top of a phase diagram taken from literature [2]. The 
compositions used span the entire Sr doping regime where the colossal magnetoresistant effect occurs. Only crystals 
were selected which have a sharp magnetic transition, and are optically flat and shiny (see the inset of Fig. S2c). 
Each cleavage surface measured with STM/S or ARPES has been investigated with low energy electron diffraction 
(LEED). Very sharp diffraction patterns were obtained of unreconstructed (lxl) tetragonal symmetry. The single 
crystallinity of the samples was confirmed by Laue measurements, of which Fig. S2d is a typical example. All of the 
above show the investigated single crystals and their cleavage surfaces to be of excellent quality [12]. 
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S 2: Sample and cleavage surface quality of the bilayer LSMO crystals, (a) Magnetization versus temperature for an x = 0.36 
sample measured using SQUID magnetometry. Data taken after zero field cooling, with an external field B = 100 G 1 1 c. The 
sample shows a sharp transition from a paramagnetic (PM) to a ferromagnetic (FM) state at Tc = 130 K with a total transition 
width of less than 5 K. (b) Magnetic transition temperatures as measured using SQUID magnetometry of LSMO with x = 0.30 
(Tjv) and 0.325 > x > 0.425 (Tc), plotted on top of the magnetic phase diagram, taken from Ref. [2]. Depicted are the onset, 
midpoint and endpoint temperatures of the transition (triangles, circles and squares), (c) LEED image of a typical bilayer 
LSMO sample, Eo = 400 eV, showing a very clear tetragonal pattern, without any signs of a 2D structural reconstruction. The 
inset shows a cleaved crystal on top of a cleavage post, with a smooth and mirror-like surface over millimeters, (d) Typical 
Laue image from our bilayer manganite crystals, indicative of the high quality and single crystallinity of the samples. 
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III. ORIGIN OF THE QP-PEAKED AREAS 
A. Fermi surface maps and T-dependent ARPES data 



For a considerable proportion of all samples investigated, ARPES energy momentum maps spanning at least one 
Brillouin zone have been measured. Upon constructing a Fermi surface by fitting peaks in the momentum distribution 
curves at the Fermi energy, /cjr-values can be extracted for each band crossing the Fermi level and compared to the 
values expected for the particular doping concentration. Representative Fermi surface maps for the QP-peaked region 
of cleaves of three different doping concentrations are shown in Fig. S3. A similar analysis can be performed on the 
non-peaked regions of the sample, constructing a 'ghost' Fermi surface by integrating the first 100 meV of the spectra. 
Figure S4a shows the doping dependence of fc values of both peaked and non-peaked regions. 

Immediately it becomes clear that the measured kp- values for the peaked regions, aside from being - in some cases 
- too many in number, do not match with those expected from band theory [13]. Moreover, there is no significant 
change in the various /c^-values as a function of doping. The highest and lowest doped compositions shown in Fig. S3 
differ 12.5% in hole doping, thus an equal change in Fermi surface area is expected between x = 0.30 and x = 0.35 
and between x = 0.35 and x= 425. From LDA calculations [13], the 2kp-values at (ir, 0) for the anti-bonding and 
bonding band are expected to change between the two extreme doping levels considered from 0.547r/a and 0.907r/a 
to 0.487r/a and 0.657r/a, respectively. The right-most panels of Fig. S3 show the fitted k positions from the FS maps 
of the three doping levels overlain on top of each other. Noticeably, the data points for the three samples fall almost 
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S 3: Large k||-area Fermi surface maps (signal integrated over Ef ± 10 meV) for bilayer LSMO with x = 0.30, 0.35 and 0.425 
taken at T = 20 K with \iv — 56 and 76 eV. The point density perpendicular to the analyzer slit is w 0.087r/a (for hv = 56). 
The large arrow in the top-left FS map indicates an unexpected electron-like Fermi surface pocket around the (n, ir) point of 
k-space. The displayed data have been interpolated with two extra slices along the angular scanning direction (i.e. along the 
(0, 0) - (it, 7r) direction). The two right-most panels show the fitted positions of the identifiable Ef crossings in the raw FS 
maps for both photon energies. Results for the various doping level are overlain to allow for easy comparison. The error bars 
in the fitted positions are smaller than the symbol size. 
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exactly on top of each other. This means that their Fermi surface areas are almost the same, despite the significant 
change in nominal hole concentration. The respective 2kp values at the (71", 0)-point are rs 0.467r/a and f=a 0.367r/a 
for the hv = 73 (nominally bonding band) and 56 eV (nominally anti-bonding band) data, respectively. Compared to 
the predictions from density functional theory, this would correspond to a doping level of x sa 0.50, thus a hole count 
some 10 to 20% higher than nominal for the doping levels discussed here. Important to note is that a doping level of 
50% corresponds to an insulating, antiferromagnetically ordered low temperature phase of LSMO, a phase which is 
not expected to support spectral weight at the Fermi level. Over many measured cleavage surfaces of different doping 
levels, some scatter (about 10%) has been observed in these respective 2k^-values of the anti-bonding and bonding 
band QP peaks, but no systematic change in the area of the Fermi surface that scales with the hole doping is evident 
in our substantial collection of high resolution ARPES data. The /c-values for the non-peaked regions on the other 
hand correspond nicely with calculations (within the uncertainty introduced by the very low spectral weight at Ef), 
and follow the expected trend as a function of doping concentration (see the blue and red symbols in Fig. S4a). 

Although often two peaks can be identified on the quasi-particle peaked areas, on multiple occasions not only single 
or double bands have been observed, but also triple and quadruple bands (see Fig. S4b). The existence of two peaks 
can readily be explained by the presence of a bonding and anti-bonding band due to the double layer of manganese 
oxide octahedra. However, it requires a lot more effort to explain more than two bands for a perfectly bilayer crystal. 
The symmetric intensity distribution of the peaks at positive and negative k^ and the switching of intensity between 
the spectra taken with different photon energy (top left panels of Fig. S3), make it highly unlikely that the multitude 
of peaks is caused by simultaneous measurement of several crystallites with different doping level or surface angles 
with respect to each other. Without having a microscopic explanation for such multi-peaked data ready and waiting, 
here we simply show such data to make the reader aware of the intrinsic variability in the spectra from the QP-pcakcd 
regions of the sample surfaces. 

In addition to comparison with calculations, the temperature dependence of the quasi-particle peaks can also give 
clues as to the origin of the QP-peaks. If the peaked spectra indeed are a true signature of the bilayer manganite, 
the presence of the peaks should be highly correlated with the transition from the insulating to poor metallic state. 
For several doping concentrations the temperature dependence of the quasi-particle weight has been examined, the 
result of which is plotted in Fig. S4c. Clearly, for none of the doping concentrations investigated do the peaks at the 
Fermi level vanish at Tq- Instead, quasi-particle spectral weight is in some cases even present to temperatures nearly 
as high as the metal to insulator transition of the cubic compound. 

Table SI shows a list of samples measured in this study. Doping concentrations spanning the entire phase space in 
which the CMR effect occurs have been measured, with multiple samples for each doping concentration. Since not 
all samples measured in this study have been spatially mapped in detail, the list is not exhaustive. At most a few 
percent of the area of each sample supports quasi-particle peaked areas in ARPES, whereas the remaining sample 
surface shows a (pseudo-)gapped spectrum, the latter being readily unified with data from bulk (optics, transport) 
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S 4: Doping and temperature dependence of the electronic structure, (a) Doping dependence of the 2fcF-values from I{E, k)- 
images recorded around the (7r,0)-point of fc-space (at T — 20 K). Round and square symbols: 2fcF~values for the inner and 
outer band of non-QP-peaked spectra. Shaded red and blue area: scatter of 2fc_F-values of the nominally bonding (hv = 73 eV) 
and anti-bonding (hv — 56 eV) bands from the QP-peaked regions of the cleaves, respectively. Yellow diamonds: 2fe_F _ values 
predicted from density functional calculations [13]. (b) I(E, k) -images at the (7r,0)-point (T — 20 K, hv = 56 and 73 eV) 
showing as many as four different QP-peaked Fermi crossings, which are clearly identified in the /(fc)-curves at Ef, shown as an 
inset to the panels, (c) Temperature dependence of the normalized QP-peaked spectral weight (integrated over Ef ± 75 meV) 
for three different doping levels. The respective bulk Tc 's of the samples are indicated. 
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Table S I: Samples measured and the estimated area-percentage supporting quasi-particle peaks per sample, as determined by 
scanning the surface in ARPES using 100 micrometer size steps. A few percent at most of the total area of a sample supports 
sharply peaked features in k-space. Only samples are shown which have been spatially mapped in detail. 



and surface sensitive spatially resolving techniques (STS) cited throughout the main text of this paper. 

B. V2xy2 in-plane symmetry for the QP-features in ARPES 

Another discrepancy between the measured Fermi surfaces from the QP-peaked regions and the predicted band 
structure for bilayer LSMO is a feature visible in all FS maps, for example that marked with a yellow arrow in Fig. S3. 
A small circular pocket can be identified, which is electron-like around (7r,7r) (and is thus markedly different from the 
band structure features discussed throughout: these all being hole-like around the (7r,7r) point). The dispersion of the 
additional electron pocket can be seen in the E(k)-image cutting through k-space along the (7r,7r) - (-7r,-7r) direction 
depicted in Fig. S5. The most intense feature in the spectrum is the main band, either the d z 2_ r 2 or the anti-bonding 
or bonding d x 2_ y 2 band that has its band bottom below BE = 1 cV at (0,0). Two other Fermi crossings, forming the 
extra electron pocket, can be seen in the image (indicated with open triangles). 

Taking a second derivative of the spectrum (in k) enhances the contrast between the main band and these features. 
Such a second derivative spectrum is depicted in the inset to the large E(k)-imagc in Fig. S5. The fitted band 
dispersions are plotted on top of the derivative spectrum and show that the extra features mirror the main band 
dispersion perfectly. Such a copy of the main band(s) can only occur if these are folded back around a smaller 
Brillouin zone boundary, resulting from a larger periodicity in real space than the tetragonal unit cell. From the 
observed dispersion of the diffraction copies, we deduce that the reconstruction of the tetragonal unit cell is most 
likely of y/2 x origin and the diffraction copy is that of the d z 2_ r 2 hole pocket originating from around (0,0), as 
illustrated in the bottom right inset to Fig. S5. 

From the above observation, it follows that QP-pcaked regions of the bilayer LSMO cleavage surface do not posses a 
lxl tetragonal unit cell. This is in sharp contrast with the bulk crystal structure of bilayer LSMO in the ferromagnetic, 
metallic low temperature phase, which is strictly tetragonal. In addition, a surface sensitive diffraction technique such 
as LEED (that is however spatially averaging over <~ 1mm 2 ) yields a strictly lxl tetragonal diffraction pattern with 
a unit cell of 3.8 A. Interestingly, no diffraction copies were ever observed on the non-QP peaked areas of a cleavage 
surface. Thus, observations of the v2 x y2 reconstructed crystal structure are confined to the QP-peaked spots only. 
Small area electron diffraction experiments in transmission electron microscopy studies on bilayered LSMO with x = 
0.40 have shown that intergrowth inclusions have a \f2 x reconstructed unit cell, with respect to the tetragonal 
structure of the bilayer matrix [15]. This would be exactly the type of reconstruction that could explain the observed 
diffraction copies seen in our high resolution ARPES data. 

C. Effect of matrix elements: photon polarization and measurement geometry in ARPES 

The ARPES and STM data presented here are argued to make a convincing case for the fact that the true spectral 
response of bilayer LSMO in the CMR region of the phase diagram is not that of a metal: pseudo-gapped, without 
QP's and thus without a Fermi surface. When arriving at such a conclusion from ARPES data it is vital to consider 
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the effects of the matrix elements, in which the choice of photon energy, polarization and experimental geometry 
can have a profound impact on the experimentally measured intensities. Already from the fact that the QP-peaked 
spectra (here reported to result from a small surface area of a typical bilayer LSMO crystal that displays stacking 
faults) have been reported to be most easily resolved using either a photon energy of hv = 56 or 73 eV [16-18], 
can one infer that matrix elements play a role for these features. This may also hold for the measurement geometry 
and polarization of the incoming photons. In fact, the ARPES studies on bilayer LSMO reported in Refs [19] and 
[20] have been taken with a measurement geometry in which the polarization vector of the incoming synchrotron 
radiation essentially lies in the sample surface plane. These two studies have shown bilayer LSMO with a doping level 
of x = 0.40 to be a 'nodal metal'. 

Therefore, we need to examine just how robust the observation (or not) of the QP-peaks is for different experimental 
conditions. Fig. S6 shows ARPES data from a QP-peaked portion of a crystal cleave with x = 0.35 (i.e. well inside 
the CMR portion of the phase diagram). In panel (a), apart from a sketch of the experimental geometry, we show 
a large area Fermi surface (FS) map, covering multiple Brillouin zones. Due to the measurement geometry adopted, 
the angle between the sample surface normal (denoted n) and the incoming light changes while recording different 
fc-space regions of the FS map. When using p-polarized light, the photon polarization vector lies for sa 89% in the 
sample surface plane for data recorded around (k x ,k y ) = (— 2n, — 2n). Moving toward a normal incidence geometry 
- i.e. for data recorded around (0, 0) - the polarization vector has equally large in- and out-of-plane contributions . 
For data taken at the (tt,tt) point of fc-space the vector has an out-of-plane contribution of 66%. On the other hand, 
for the same geometry, the polarization vector is 100% in-plane, irrespective of the probed area of fc-space, for data 
recorded with s-polarized light. Hence, the data contained in the FS map shown in Fig. S6a, taken together with 
s-polarization, offer an ideal test of the robustness of the presence of QP-peaked against polarization and geometry 
induced matrix element effects. 

Panel b of Fig. S6 shows l(E, k) images that cut fc-space around (— tt, —tt) and (— 2-7T, — 2tt), taken with both s- and 
p-polarized light. The amount of in-plane polarization varies between 100 and 66%, as marked in the figure. Despite 
this, in all four cases the bands comprising the X-centered hole-like FS, here cut along the zone diagonal, clearly 




S 5: The QP-peaked areas of a cleavage surface show a y/2xy/2 reconstruction. I(E, k) image obtained from an x = 0.375 
sample (hv = 56 eV, T = 20 K) showing a copy (indicated with open triangles) of the main band (indicated with a closed 
triangle). The I(E,k) image was recorded along the (0, 0)— (7T, tt) direction in k-space, as indicated with a thick purple line in 
the small FS map (top right). The inset to the I(E, k) image shows a second derivative of a similar spectrum for x=0.30, with 
the red lines showing the results of a fit to the peaks in the I(k) curves for the main band and one branch of the diffraction 
replica. The bottom right schematic FS illustrates how folding of the bands due to the presence of a x/^xv 7 ^ reconstructed 
(smaller) Brillouin zone can lead to the observed band with its band bottom at (n, tt). 



7 




S 6: QP-peaked ARPES data from bilayer LSMO (a; = 0.35, T — 20 K) recorded using different measurement geometries and 
photon polarizations, (a) Top: sketch of the measurement geometry, whereby 'n' is the surface normal and 'e' the outgoing 
photoelectron beam on its way to the analyzer slit. The black double-headed arrows indicate the polarization of the synchrotron 
radiation: either perpendicular (s) or parallel (p) to the plane spanned by the incoming light, the outgoing photoelectrons and 
the sample normal. Panel (a) bottom: large area Fermi surface map of bilayer LSMO (p-polarized light) recorded by rotation 
of the sample around the axis indicated in the sketch. A schematic (plan view) of the measurement geometry is shown for data 
taken through the (0, 0) fc-space location. Note that the in-plane component of the light polarization increases upon rotation 
of the sample so as to reach the (— 2tv, —2n) point. The cuts made through the dataset that are shown in panels (b) and (c) are 
indicated using colour-coded triangles, (b) 1(E, fc)-images from two different Brillouin zone diagonal cuts in fc-space through 
the main hole-like bands, measured both with p- and s-polarization, including sketches of the measurement geometry. The 
in-plane component of the polarization of the incoming photons is marked on each I(_E, k) image as a percentage, (c) Three 
I(-E, fc)-images from fc-space cuts through the 'zone face' region of the Fermi surface, recorded with p-polarization. Again, 
the in-plane polarization component of the photon polarization is given, as well as a sketch of the measurement geometry for 
sub-panel c3. The right-most sub-panel shows I(i?)-curves for the three 1(E, fc)-images in panel c and the one in panel b4, 
obtained by integration of a small fc-range around &f (as indicated with the red dotted box in sub-panel c3) and offset for 
clarity. 

show QP-peaks at Ep (indicated by colored triangles whose colour matches the symbols marking the k-space cuts in 
the Fermi surface map). 
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Additionally, the cut shown in panel bl shows two electron-like bands centered around (— 2tt, —2tt), one of which 
is very shallow and has its band bottom below a binding energy of BE = 50 meV. Panel b3 shows a copy of these 
features centered around (—n, —n), due to the decreased in-plane crystal symmetry of the stacking fault intergrowths, 
as has been discussed in connection with Fig. S5. 

The same cut in fc-space taken with s-polarized light (panel b2) shows a strong suppression of the two electron-like 
features due to the selection rules. This also true for the (— n, — 7r)-located diffraction copy of these features in panel 
b4. 

The 'extinction' of the two T-centered electron pockets (and their diffraction copies) upon use of s-polarized light 
indicates that they have mainly d z i_ r i character, in keeping with predictions from LDA calculations [13, 14] . Due to 
the measurement geometry, the selection rules only apply strictly at center of the analyzer slit (i.e. the centre of the 
k-space axis of each l(E, k) image), as only this point lies in the mirror plane defined by the incoming photons, the 
surface normal and the outgoing photoelectrons. This hinders an unambiguous determination of the orbital character 
of the bands comprising the A-centered hole pocket. 

Irrespective of the selection rules, the QP-peak spectral weight of the hole-like features is almost constant compared 
to the incoherent spectral weight at higher binding energies across all four panels of Fig. S6b. In particular, the 
similarity between panels b3 and b4 shows that the zone-diagonal QP-peak is still robustly observed, despite a 
significantly difference in in-plane polarization for the two 'nodal' cuts. 

This observation also holds for the QP-peaks observed around the zone face regions of /c-space, see Fig. S6c. In 
this panel, three I(E, fc)-images are presented that have a decreasing in-plane contribution of the photon polarization. 
Cut cl, with the highest in-plane contribution, docs show reduced QP spectral weight compared to the other two 
cuts. However, it still has significant spectral weight at Ep and is thus far from being 'pseudo-gapped'. This can also 
be seen in the I(_E)-curves for k ~ kp depicted in the rightmost panel of Fig. S6c, which also includes a zone-diagonal 
k F I(.E)-trace for comparison. 

This section shows that although the measurement geometry does play a role in the intensity of the QP-peaks at 
Ep, the observation of QP-peaks at the various /c-locations on the Fermi surface is not restricted to a particular set of 
experimental conditions, but is restricted to particular - minority - locations on the sample surface from which signals 
from intergrowth patches are picked up in the experiment. Also important is the fact that the lack of QP-peaks 
from the majority of the cleavage surface is not a particular result, but generally applicable, also for experimental 
geometries with a large in-plane component of the photon polarization, and for a wide range of photon energies. The 
modest sensitivity of the data from both the non-QP-peaked and QP-peaked regions to experimental conditions only 
serves to strengthen our conclusion that the general spectral signature of bilayer LSMO in the ferromagnetic region 
of the phase diagram is that of a Fermi-surface-less non-metal, without coherent spectral weight at Ep, and that the 
QP-features and Fermi surfaces we observe originate from N^2 intergrowths. 
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IV. PSEUDO-PERIODIC TEXTURES SEEN WITH STM 



The pseudo-periodic texture as seen with STM on the bilayer manganites has been observed on samples with 
various doping concentrations. Fig. S7 shows topographies for four different samples supporting such textures. Auto- 
correlation traces of these four topographs are shown in Fig. S7e, in all cases having correlation peaks at distances 
ranging from ~ 5-10 unit cells. There is no simple trend in these characteristic distances as a function of doping, and 
nor are the periodic structures pinned to the real-space lattice vectors. 




S 7: (a)-(d) Pseudo-periodic texture observed in the topographic STM images from different Sr doping concentrations of bilayer 
LSMO. The setup conditions are similar to those in Fig. 1 of the main paper. The vertical corrugation of these semi-regular, 
long wavelength features is, at less than 1 A too small to be due to differing atomic layers of the crystal. Doping levels are 
given in the legend to panel (e), which shows line traces through the autocorrelation plots derived from the topographs. These 
traces show that the data support differing correlations of range between 5-10 unit cells. Each autocorrelation figure is shown 
as an inset to the corresponding STM topograph. 
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